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Abstract. The orbits and physical parameters of three close, double-lined GO eclipsing binaries have been derived combining 
H P , V T ,B T photometry from the Hipparcos/Tycho mission with 8480-8740 A ground-based spectroscopy. The setup is mim- 
icking the photometric and spectroscopic observations that should be obtained by GAIA. The binaries considered here are all 
of GO spectral type, but each with its own complications: V781 Tau is an overcontact system with components of unequal tem- 
perature, UV Leo shows occasional surface spots and GK Dra contains a S Scuti variable. Such peculiarities will be common 
among binaries to be discovered by GAIA. We find that the values of masses, radii and temperatures for such stars can be 
derived with a 1-2% accuracy using the adopted GAIA-like observing mode. 

Key words. surveys:GAIA - stars:fundamental parameters - binaries:eclipsing - binaries:spectroscopic 



1. Introduction 

GAIA is a challenging Cornerstone mission re-approved by 
ESA last May for a lunch by around 2010. It is aimed to provide 
micro-arcsec astrometry, 10-band photometry and medium res- 
olution 8480-8740A spectroscopy for a huge number of stars, 
with completeness limits for astrometry and photometry set 
to V — 20 mag. Each target star will be measured around 
a hundred times during the five year mission life-time, in a 
fashion similar to the highly successful operational mode of 
Hipparcos. The astrophysical and technical guidelines of the 
mission are described in the ESA's Concept and Technology 
Study (ESA SP-2000-4), in the papers by Gilmore et al. (1998) 
and Perryman et al. (2001), and in the proceedings of con- 
ferences devoted to GAIA and edited by Straizys (1999), 
Bienayme and Turon (2002), Vansevicius et al. (2002) and 
Munari (2003). 

In Paper I of this series, Munari et al. (2001), we have 
started to provide reasonable orbits for a number of new eclips- 
ing binaries and to evaluate expected performances of GAIA on 
eclipsing binaries with an emphasis on the achievable accuracy 
of derived fundamental stellar parameters like masses and radii. 
The expected number of eclipsing binaries to be discovered by 
GAIA is ~ 4 x 10 5 . Some 10 5 of these will be characterized 
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as double-lined in GAIA spectral observations. This is a huge 
number, many orders of magnitude larger than the total of SB2 
eclipsing binaries so far investigated from ground-based ob- 
servations (cf. Andersen 1991; Batten, Fletcher & MacCarthy 
1989). Perhaps the orbits and stellar parameters could be de- 
rived from GAIA observations at a few percent error only for 
a few percent of them. But this still represents a two-orders 
of magnitude increase compared to all ground-based observ- 
ing campaigns during the last century. Data obtained by GAIA 
should be able to provide reasonable solutions as ground-based 
follow-up campaigns will be very time consuming. It is there- 
fore of great interest to investigate the expected performances 
of GAIA on eclipsing binaries. The purpose of this series of pa- 
pers is to contribute to the fine tuning of the last details in the 
mission planning as well as to define the strategy to analyze 
the massive spectroscopic and photometric data flow on eclips- 
ing binaries that is completely unprecedented. In the meantime, 
this series of papers will focus on eclipsing binaries unknown 
or poorly studied in the literature so far. 

Paper I outlines the framework of the project and adopted 
methods, and the reader is referred to it (and the references 
therein) for details. In short, Hipparcos/Tycho photometry is 
adopted as a fair simulation of typical GAIA photometric 
data. The satellite spectroscopic data is simulated by devoted 
ground-based observations obtained with the Asiago 1.82m + 
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Table 1. Programme eclipsing binaries. Data from the Hipparcos Catalogue. H P , B T , V T are out-of-eclipse median values. 

Name Spct. H P B T V T a J20 oo S J2 ooo parallax dist /i a fi 6 

(hms) (° ' ") (mas) (pc) (masyr -1 ) (masyr -1 ) 

V781 Tau HIP 27562 GO 8.71 9.41 8.74 05 50 13.12 +26 57 43.4 12.31 + 1.35 81*> -0.084 + 0.004 -0.091 + 0.004 

UVLeo HIP 52066 GO 9.20 9.78 9.00 10 38 20.77 +14 16 03.7 10.85 + 1.16 91'° 3 -0.007 + 0.004 0.010 + 0.004 

GKDra HIP 82056 GO 8.92 9.19 8.81 16 45 41.19 +68 15 30.9 3.37 + 0.69 297™ -0.013 + 0.002 0.014 + 0.002 



Echelle + CCD, set up to mimic the expected GAIA spectra. 
Precision of the results of our investigation can be considered 
as a lower limit to the accuracy obtainable from GAIA at the 
given source S/N, because (a) GAIA will observe in many more 
photometric bands than Hipparcos/Tycho and with far higher 
accuracy even in the narrow bands, thus both increasing light- 
curve mapping as well as accuracy of information on stellar 
temperature, limb-darkening and reddening; and (b) GAIA will 
acquire at least twice as many spectra per star than considered 
here due to obvious limitations in the telescope access time. 

2. Target selection 

Similar to Paper I we have selected both some brand-new 
eclipsing binaries (i.e. without a spectroscopic or photometric 
orbit solution in the literature) as well as binaries with already 
published orbital solutions (however not in the GAIA spectral 
range) that can serve as an external comparison. Their basic 
properties are quoted in Table 1 . 

V781 Tau. This is a GO over-contact (~23%) binary (P ~ 
0.4 days) with stars of unequal temperature. It is known to un- 
dergo period changes (Donato et al. 2003, in preparation), inter- 
preted by Liu and Yang (2000) as shrinkage of the secondary. 
A spectrophotometric orbit of moderate quality has been pub- 
lished by Lu(1993). 

UV Leo. This is a GO short period binary (P = 0.6 days) 
showing intrinsic variations caused by cool spots on the sec- 
ondary component (cf. Mikuz et al. 2002). Orbital parameters 
have been derived from UBV photometric data by Frederik & 
Etzel ( 1 996) and from 4430-6800 A spectroscopic observations 
by Popper (1997). 

GK Dra. This is a newly discovered eclipsing binary, the 
only existing information in the literature being BV photomet- 
ric monitoring by Dallaporta et al. (2002). The authors showed 
that the photometric period listed in the Hipparcos Catalogue 
(~ 17 days) is wrong (the actual one being 9.97 days), and that 
the secondary star has intrinsic variability of a 5-Sct type. 

3. Observations 

As explained above we use Hipparcos photometry as a lower 
limit to the photometric information expected from GAIA. The 
accuracy of Hipparcos photometry is lower, but the number of 
observations of each star with only a limited number of points 
sampling the eclipses is similar. Table 2 gives details on the 
number of observations of each star and their accuracy. 



All spectral observations were obtained in the same mode 
as in Paper I, i.e. at 0.25 A/pix dispersion and ~0.50 A res- 
olution over the 8480-8740 A wavelength range (therefore a 
resolving power R = A/ A A = 17,000). 

The spectroscopic observations have been collected with 
the Echelle+CCD spectrograph on the 1.82 m telescope op- 
erated by Osservatorio Astronomico di Padova atop of Mt. 
Ekar (Asiago). A 2.2 arcsec slit width was adopted to meet 
the resolution requirement. The detector has been a UV coated 
Thompson CCD with 1024x1024 square pixels of 19jUm size. 
The GAIA spectral range is covered without gaps in a single 
order by the Asiago Echelle spectrograph. The actual observa- 
tions however extended over a much larger wavelength interval 
(4550-9600 A). Here we will limit the analysis to the GAIA 
spectral interval; the remaining, much larger wavelength do- 
main will be analyzed elsewhere together with devoted multi- 
band photometry from ground based observations. The spectra 
have been extracted and calibrated in a standard fashion using 
IRAF software packages running on a PC under the Linux op- 
erating system. The high stability of the wavelength scale of the 
Asiago Echelle spectrograph has been discussed in Paper I. The 
results of radial velocity measurements are given in Table 3. 

4. Modeling 

We use an upgrade of the setup described in Paper I. The bi- 
nary modeling code (Wilson 1998) was combined with Van 
Hamme's limb darkening coefficients (van Hamme 1993), a 
fitting package, a graphical user interface and utilities like red- 
dening corrections to form PHOEBE (Prsa 2003). The pack- 
age is able to run on any Unix platform. It may constitute the 
first step toward automated solution-finding routines that will 
be needed to interpret the vast number of binary systems to be 

Table 2. Number of Hipparcos (H P ) and Tycho (B T , V T ) photo- 
metric data and ground based radial velocity observations, their 
mean S/N and standard error for the three programme stars. 
Error for radial velocity is in km s - ' . 

Hip Tyc RV 



N g-(Hp) N cr(B T ) o-(W) N S/N cr(RV) 

V781 Tau 61 0.014 81 0.18 0.15 41 35 8 

UVLeo 96 0.015 150 0.21 0.17 29 30 10 

GKDra 124 0.017 179 0.15 0.15 35 45 3 



T. Zwitter et al.: Evaluating GAIA performances on eclipsing binaries. 



3 



Table 3. Journal of radial velocity data. The columns give the spectrum number (as from the Asiago 1.82 m Echelle+CCD log 
book), the heliocentric JD, and the heliocentric radial velocities (in km s _1 ) for both components. An asterisk marks the spectra 
with a too severe blending of the lines for a meaningful measurement of radial velocities of each component. The latter have not 
been used in modeling of the binaries. 
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observed by GAIA. All results were independently derived also 
by the WD98K93 code (Milone et al. 1992) and WD2002 code 
(Kallrath et al. 1998) that are briefly described in Paper I. We 
found that the results are in agreement. 

The usual approach to binary star modeling is to use only 
relative photometry obtained in each filter. Depths of eclipses 
in different filters constrain the ratio of the stellar temperatures, 
while the absolute temperature scale is tuned by judging the 
primary star temperature from the system colour. 

In our case both stars are of similar brightness and the light 
curves are quite noisy. This requires some modifications to the 



usual approach. Hipparcos observed in three filters. The obser- 
vations obtained in the broad band H P filter have an accept- 
able accuracy, while those in the Tycho experiment's Bj and 
Vt bands are generally very noisy. We use the absolute system 
colours at quarter phase to fix the absolute temperature scale. 
The transformation between the Tycho and Johnson systems is 
the same as in Paper I: 

Vj = V T - 0.090 x (B - V) T (1) 
(B - V)j = 0.85 x (B - V) T (2) 



4 



T. Zwitter et al.: Evaluating GAIA performances on eclipsing binaries. 



Table 4. Modeling solutions. The uncertainties are formal mean standard errors to the solution. The last four rows give the 
standard deviation of the observed points from the derived orbital solution. 



parameter(units) 


V781 Tau 




UV Leo 




GK Dra 




Period (days) 


0.34490857 


+ 


0.0000001 


0.600086 ± 


0.000001 


9.9742 


+ 


0.0002 


Epoch (HJD) 


2447962.46572 


+ 


0.00016 


2448500.560 ± 


0.001 


2452005.56 


+ 


0.03 




2.4478 


+ 


0.002 


3.957 ± 


0.087 


28.92 


+ 


0.35 


V y (km s _1 ) 
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0.10 


3.9 ± 
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0.24 


Ri (Ro) 


0.759 




0.007 


0.973 ± 


0.024 


2.431 




0.042 




1.111 




0.007 


1.216 ± 


0.043 


2.830 




0.054 


ivi i \iyiQ) 


0.510 


+ 


0.006 


1.210 + 


0.097 


1.460 


+ 


0.066 




1.150 


+ 


0.027 


1.110 + 


0.100 


1.810 


+ 


0.109 


Mj, ; i 


4.950 


+ 


0.025 


4.590 ± 


0.094 


1.960 


+ 


0.075 




4.280 


+ 


0.020 


4.390 ± 


0.113 


1.770 


+ 


0.072 


loggi (cgs) 


4.380 


± 


0.012 


4.540 ± 


0.053 


3.830 


+ 


0.033 


log ft (cgs) 


4.410 


+ 


0.016 


4.310 ± 


0.055 


3.790 


+ 


0.041 


o-Rvx2 (km s _1 ) 


13.8 






17.6 




2.71 






cr(5 ( ) (mag) 


0.193 






0.228 




0.187 






cr(V t ) (mag) 


0.173 






0.227 




0.199 






o-(ff p ) (mag) 


0.020 






0.028 




0.028 







Temperatures of both stars are similar, so the temperatures of 
the stars, T\, T2 and their radii Ri, R2 are connected to the 
surface-weighted effective temperature of the source at quar- 
ter phase T 1+ 2 by the relation: 

R 2 J 4 1+ R 2 2 T* = (R]+R 2 2 )T 4 l+2 (3) 

First the Tycho colour index at quarter phase of the model fits 
to the Bj and Vj light curves was transformed to the Johnson 
system (Eq. (2)) and the effective temperature r [+ 2 was deter- 
mined. Modeling of the better quality H P band observations 
yielded the temperature ratio and, by use of Eq. (3), also the 
absolute temperatures of the two stars. The process was reiter- 
ated several times to reach a self-consistent solution. 

Some colour calibrations proposed recently (Bessell 2000) 
differ from Eq. (2) and cause effective temperature offsets of 
~ 100 K. We will comment on the changes of the results if 
these relations were used in the Discussion. 

5. Results 

Table 4 quotes the derived system parameters together with 
their formal errors. Table 5 compares the derived distances to 
the astrometric results from Hipparcos. The data and the curves 
from the model solutions are plotted in Figs. 1-3. 

We note that model fits are generally acceptable. The dif- 
ferences are chiefly due to noise in the data and to some de- 
gree due to intrinsic variability of the stars. A limited number 



of epochs and their long timespan make modeling of transient 
phenomena such as stellar spots unfeasible. This will generally 
be also the case with data obtained by GAIA. The results were 
obtained assuming the stars are co-rotating. Next we discuss in 
turn the results for each of the objects. 

5.1. V781 Tau 

V78 1 Tau is an overcontact binary with different primary and 
secondary temperatures. Light curve modeling fixes the quarter 
phase magnitudes (O = 0.75) to V T = 8.56 and B T = 9.16. 
This corresponds to the colour index (B - V)t = 0.60 or 
(B - V)j = 0.51 (Eq. (2)) which gives T i+2 = 6240 K. 
This result was used to constrain the temperatures of the two 
stars through Eq. (3). Note that the magnitudes quoted in the 
Hipparcos catalogue (Table 1) would give somewhat different 
colours. However these magnitudes are just a suitable mean of 
all observations, also the ones close to the photometric eclipses. 
Therefore it is correct to use the quarter phase light curve fit and 
not the mean colours. 

Spectroscopic observations determine absolute size of the 
system and individual masses as a function of the system incli- 
nation. A detailed reflection treatment was used to compute the 
photometric curves. The H P light curve constrains relative sizes 
and temperature ratio of both stars. We found the system is ac- 
tually filling its Roche lobes up to the L\ point. The stars are of 
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Fig. 1. Hipparcos Hp and Tycho W, Br, (B - V)t lightcurves 
of V781 Tau folded onto the period P = 0.34490857 days. 
Radial velocity measurements in the GAIA spectral interval 
from Table 3 are given on the top, with '+' signs marking 
blended spectral lines that were not used for modeling. The 
curves represent the solution given in Table 4. 
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Fig. 2. Hipparcos and Tycho Vt,Bt,(B-V)t lightcurves of 
UV Leo folded onto the period P = 0.600086 days. Radial ve- 
locity measurements in the GAIA spectral interval from Table 3 
are given on the top, with '+' signs marking measurements 
around primary eclipse that were not considered in modeling. 
The curves represent the solution given in Table 4. 



unequal temperature (T1-T2 ~ 220 K). This difference was ex- 
plained by mass transfer between the stars and the correspond- 
ing gravitational energy release (Liu & Yang 2000). A small 
period decrease (dP/P = -5.0 x 10~ n ) was also claimed to be 
an effect of mass transfer. We note that any mass lost from the 
system through the L 2 point would carry away roughly twice 
the mean value of the specific angular momentum. Mass loss 
through the L2 point can therefore decrease the total angular 
momentum of the system, so it may be partially responsible for 
shortening of the orbital period. The value of the time deriva- 
tive is too small to be detectable from data used in this study. 



Lu (1993) published a spectrophotometric study roughly at 
the same accuracy level as reported in Table 4. The values of 
individual parameters are generally consistent, with some dif- 
ferences possibly arising from the simplified software he used 
for modeling. In particular he adopted lower effective tempera- 
tures (T h2 =5950, 5861 K) but with a large error bar of 200 K. 
Therefore the system in his analysis turns out to be fainter and 
at a smaller distance (72 pc). 

We note that the formal error bars on temperatures as given 
by the WD98 code can be increased due to systematic effects. 
True uncertainty can reach 100 K, increasing the uncertainty on 
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Fig. 3. Hipparcos Hp and Tycho Vt, Bj, (B - V)p lightcurves of 
GK Dra folded onto the period P = 9.9742 days. Radial veloc- 
ity measurements in the GAIA spectral interval from Table 3 
are given on the top. The curves represent the solution given in 
Table 4. 



the distance (Table 5) to 4.5% or 3.6 pc. Temperatures of both 
stars may be also influenced by reddening. V 781 Tau lies just 
0.2° from the Galactic plane. One may expect E(B - V) = 0.09, 
and Ay = 0.3 mag (Perry and Johnston 1982). In our calibra- 
tion the effective temperature T\ + 2 would raise to 6540 K and 
the bolometric magnitude of the system would be brighter by 
0.33 mag. Note that this brightening almost cancels out with 
the value of the total extinction. So reddening has little influ- 
ence on the distance of the system reported in Table 5. 



5.2. UV Leo 

UV Leo is a close binary with a pronounced spot activity that 
is expected to be common between G/K type binaries to be ob- 
served by GAIA. The spots cause vertical offsets in the bright- 
ness of the object on a time-scale of weeks to months (Mikuz 
et al. 2002). Such intrinsic variability may be contributing to 
the scatter of Hp observations in Fig. 2. Magnetic activity may 
be also responsible for part of the scatter of the radial veloc- 
ity curves (<Xpy 1>2 = 17.6 km s _1 , Table 4). In fact the Ca II 
lines from the secondary on JD 2451896 show hints of multi- 
component profiles, typical for spotted stars. This structure, 
though below the level suitable for detailed analysis in our (and 
usually also GAIA's) coverage of the Ca II lines, obviously in- 
creases the scatter of derived radial velocities. 

The fits to the Vp and Bp curves give a quarter phase colour 
(B - V) T = 0.72, corresponding to (B - V)j = 0.61. This is 
consistent with the colours derived by Popper (1997). For main 
sequence stars this colour index translates into T l+ 2 = 5900 K. 
This constraint was adopted during our spectrophotometric 
model fitting. 

Popper (1997) published a spectrophotometric solution de- 
riving the average masses, radii and temperatures of both stars. 
Here we derive the parameters also for individual stars. The 
results are generally consistent. 

5.3. GK Dra 

Similar to UV Leo, GK Dra also features intrinsic variability 
of its components. The variability is however not caused by 
spots but by a likely 5-Sct variability on the secondary star 
(Dallaporta et al. 2002). This variability has an amplitude of 
~ 0.05 mag, so it is partially responsible for the scatter in the 
Hp curve in Fig. 3. The Vp and Bp curves are very noisy. Still 
they provide an average quarter phase colour (B - V)p = 0.39, 
corresponding to (B - V)j = 0.33 and effective temperature 
Ti+2 = 7000 K. The photometry to be obtained by GAIA will 
be of much higher accuracy (cr ~ 0.001 mag) than Tycho ob- 
servations. This will provide for accurate colour information 
also during eclipses and therefore constrain the temperature of 
either star. 

Hipparcos catalogue lists an orbital period of 16.96 days. 
Dallaporta et al. (2002) showed by a devoted ground-based ob- 
servation campaign that the true period is 9.97 days. The error 
in the Hipparcos results can be traced to the fact that the orbital 
period had to be derived from only 124 points. The system is 
detached so only 15 point fell into either eclipse. Spectroscopic 
information obtained by GAIA will greatly alleviate such prob- 
lems (see Zwitter 2003 for detailed simulations). This is a con- 
sequence of the fact that every radial velocity point contributes 
to period determination and not only those falling into eclipses 
as for photometric observations. 



6. Discussion 

Our analysis used the Tycho to Johnson colour transformation 
from the Hipparcos catalogue as given in Eqs. (1) and (2). The 
magnitude measurements themselves were obtained from the 
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Hipparcos and Tycho-1 epoch photometry as available through 
the CDS. Recently Bessell (2000) published modified calibra- 
tions that would make the (B - V)j colours redder by 0.03 to 
0.04 mag. The T 1+2 temperatures for V781 Tau, UV Leo and 
GK Dra would be lower for 120 K, 100 K and 160 K, respec- 
tively. A modified Vj-Vt vs. (B-V)t relation would also make 
their apparent Vj magnitudes ~ 0.01 mag brighter. Such small 
corrections cannot significantly modify the limb-darkening and 
other coefficients that depend on the absolute value of the tem- 
perature. But they do change the bolometric magnitudes and 
so distances. In our case the absolute bolometric magnitudes 
for V781 Tau, UV Leo and GK Dra would be 0.28, 0.19 and 
0.46 mag fainter and the derived distances 14, 9 and 23 % 
larger. The issue of absolute colour calibrations of the Tycho 
passbands does not seem to be a closed one. The new version of 
the Tycho catalogue (Tycho-2) quotes the old calibration (Eq. 
(2)) again. We therefore prefer to remain with the same calibra- 
tion as used in Paper I with the possible modifications clearly 
spelled out. 

7. Conclusions 

The paper clearly demonstrates the potential of GAIA to derive 
accurate orbital solutions even for stars with intrinsic variabil- 
ity or for contact cases. GAIA will observe any object only 
around a hundred times. This will complicate the determina- 
tion of orbital period of wide detached systems. Spectroscopic 
information will be particularly useful to determine the or- 
bital period in such cases and also for a vast majority of bi- 
naries which are non-eclipsing. Spectroscopic information can 
be used also to derive orbital eccentricity as demonstrated by 
GK Dra. 

Absolute scale of the system provided by spectroscopic or- 
bit can be used to derive masses and sizes of the system compo- 
nents at a 1-2% level (Table 4). So these stars can be absolutely 
placed on an H-R diagram. Exact coevality of both stars in a 
binary make for a useful study of stellar isochrones. Munari 
(2003) discusses how additional information, like metallicity, 
will be obtained from the GAIA data. 

The distances derived from orbital solutions compete or are 
superior to the Hipparcos astrometric measurements. We note 
that the present analysis may be influenced by uncertain cali- 
brations in the noisy photometry obtained from the Hipparcos 
Tycho experiment. But for the case of GAIA the errors quoted 
in Table 5 are realistic, as the stellar temperatures and redden- 
ing will be known with high precision from a multi-band pho- 
tometry. Note also that measurement of distances from orbital 
solutions, especially for overcontact binaries, is limited only by 
relative faintness of the objects at large distances. So hot con- 
tact binaries will be a useful tool to gauge distance throughout 
the Galaxy and beyond. 

GAIA will be able to detect also intrinsic variability of bi- 
nary components. Degree of derivable physical information de- 
pends on the nature of the variability. Stellar spots will be very 
common but difficult to describe. These are transient phenom- 
ena, so the star will look different on each of the 100 transits 
during the 5-yr mission lifetime. This can be seen also in our 
data. Different levels of quarter phase maxima in the V78 1 Tau 



Table 5. Comparison between the Hipparcos distances and 
those derived from the parameters of the modeling solution in 
Table 4. Only formal errors quoted in Table 4 were taken into 
account. As explained in the text the actual distances may have 
a bit larger uncertainties. 

Hipparcos this paper 
(pc) (pc) 

V781 Tau 81^ 81+1.0 

UV Leo 91 >f 92+6 

GK Dra 297™ 313+14 



light curve (Cereda et al. 1988) were used to claim the presence 
of polar spots (Lu 1993). But Hipparcos light curves do not re- 
veal such details. Also UV Leo is an object with occasional 
spots that change the overall system brightness. The fact that 
we ignored such phenomena but still derived quite accurate or- 
bital solutions in two systems suggests that magnetic phenom- 
ena cannot jeopardize the derivation of binary star parameters 
to some limit of accuracy. Other types of variability, like (5-Sct 
variability in GK Dra (Dallaporta et al. 2002, Zwitter 2003) 
maintain its phase, so they will be easily detectable from GAIA 
data. Orbital period changes, e.g. due to passages of the third 
body will be quite uncommon and difficult to detect due to a 
limited mission lifetime. 

This work reassures us of the high quality of physical in- 
formation recoverable from GAIA's observations of eclipsing 
binaries. In future papers of this series we plan to explore more 
objects with intrinsic variability as well as some double lined 
systems with triple components. 
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